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Abstract 

We propose a method to solve the Non Perturbative Renormalization Group equations for the n- 
point functions. In leading order, it consists in solving the equations obtained by closing the infinite 
hierarchy of equations for the n-point functions. This is achieved: i) by exploiting the decoupling 
of modes and the analyticity of the n-point functions at small momenta: this allows us to neglect 
some momentum dependence of the vertices entering the flow equations; ii) by relating vertices at 
zero momenta to derivatives of lower order vertices with respect to a constant background field. 
Although the approximation is not controlled by a small parameter, its accuracy can be systemat¬ 
ically improved. When it is applied to the 0{N) model, its leading order is exact in the large N 
limit; in this case, one recovers known results in a simple and direct way, i.e., without introducing 
an auxiliary field. 
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1 Introduction 


In the last ten years a considerable amount of effort has been devoted to the study of 
the Non Perturbative Renormalization Group (NPRG), and to the development of new 
approximation schemes to solve the corresponding inhnite hierarchy of equations for the n- 
point functions [1-4], In this context, the so-called “derivative expansion” has been widely 
used. It dehnes a systematic approach that has led to many successful applications in a 
variety of domains [5-7]. However, the derivative expansion is a good approximation to 
n-point functions only when the external momenta are smaller than the lowest mass in the 
problem. In particular, for massless theories, the derivative expansion provides information 
only on the n-point functions (and their derivatives) at zero momenta. This is not enough 
for applications which require the knowledge of the full momentum dependence of the 
n-point functions. In these cases, a new approximation scheme is necessary. 

To our knowledge, all efforts in this direction [9] have been based on various forms 
of the early proposal by Weinberg [8]: one truncates the inhnite tower of how equations 
for the n-point functions considering only vertices up to a given number of legs, possibly 
using various ansatzs for some of them. This leads to approximations similar to those 
used when solving Schwinger-Dyson equations [10]. However, despite the fact that very 
encouraging results have been obtained problems remain with truncations around zero 
external helds [11]. 

It is then interesting to look for alternative schemes to solve NPRG at hnite external 
momenta and it is the purpose of this Letter to present a method to calculate correlation 
functions at arbitrary momenta within the NPRG. As it is the case for the derivative 
expansion, the proposed strategy takes into account, at each order, an inhnite number of 
vertices. The accuracy is expected to be comparable to that achieved with the derivative 
expansion when it is used to calculate the ehective potential. The method exploits two 
properties of the NPRG: the decoupling of short wave-length modes, and the analyticity 
of the n-point functions at small momentum (guaranteed by the infrared regulator), in 
order to neglect some of the momentum dependence in the vertices which enter the how 
equations. Then, it becomes possible to close the inhnite hierarchy of equations by cal¬ 
culating vertices at zero momenta as derivatives of lower order vertices with respect to 
a constant background held. Although the approximation is not controlled by any small 
parameter, its accuracy can be systematically improved, in a way similar to what can be 
done in the derivative expansion. 

The leading order of the proposed approximation scheme shares some of the nice 
properties of the leading order of the derivative expansion, the so-called local potential 
approximation (LPA): it is exact at one loop order, and it also reproduces the exact large 
N limit of the 0{N) model. Of course, the LPA provides only information on momentum 
independent quantities, such as the ehective potential, whereas the present method does 
not have this limitation. In fact, as an illustration, we use it to solve the NPRG equations 
for the full momentum dependent n-point functions of the 0{N) model in the large N 
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limit. This can be done directly, i.e., without having to introduce an auxiliary held. To 
our knowledge, this is the hrst time that, within the NPRG, such a calculation is done 
for this model. 

The method presented here is an improvement of the one that we have applied in a 
previous paper to the calculation of the transition temperature of a weakly repulsive Bose 
gas [12]. Compared to the method proposed in [12], the present one is conceptually simpler 
as it involves a single approximation. Its numerical implementation will be discussed in a 
separate publication [13]. 


2 The NPRG and the derivative expansion 


Let us start by recalling some basic features of the NPRG. Although most of the 
arguments in this paper have a wider range of applicability, we shall consider a scalar 
held theory dehned by an Euclidean action R in d dimensions, of the generic form: 

S = J + V{(p)^ , ( 1 ) 

where V{ip) is a polynomial in ip. It is understood that the parameters of the action 
(hidden in V{ip)) and the held normalization are hxed at a microscopic momentum scale 
A (which may be inhnite). The NPRG equations relate the classical action to the full 
ehective action. This relation is obtained by controlling the magnitude of long wavelength 
held huctuations with the help of an infrared cut-oh, which is implemented [1-4] by adding 
to the classical action (1) a regulator of the form 

= ^j (2) 

where denotes a family of “cut-oh functions” depending on a parameter k. The role 
of AS'k is to suppress the huctuations with momenta p ^ n, while leaving unahected the 
modes with p^ n. Thus, typically Rk{p) —> when p k,, and Rnip) —^ 0 when p^ k. 
There is a large freedom in the choice of Rk{p), abundantly discussed in the literature 
[14-17]. 


We denote the ehective action in the presence of the regulator by rK[0], where 0 is the 
average held, (j){x) = {<p{x)). When k, ^ A quantum huctuations are suppressed and Pa[0] 
coincides with the classical action. As k decreases, more and more quantum huctuations 
are taken into account and, as k —0, Pk=o[0] becomes the usual ehective action r[0]. 
In other words, as k, decreases from A to 0, P^ interpolates between the classical action 
and the full ehective action (see e.g. [5]). The variation with k, of rK[0] is governed by the 
following how equation [1-4]: 
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Fig. 1. The two diagrams contributing to the flow of the 2-point vertex. The lines represent 
dressed propagators, G^. The cross represents an insertion of df^Rk- The vertices denoted by 
black dots are T^f^ and 

where is the second derivative of w.r.t. 0 (see Eq. (5) below). Eq. 3 is the master 
equation of the NPRG. Its right hand side has the structure of a one loop integral, with 
one insertion of 

As well known [18], the effective action r[0] is the generating functional of the one- 
particle irreducible n-point functions. Similarly, for a given value of k, we dehne the 
n-point functions F^^h 

/ r 

(4) 

By deriving eq. (3) with respect to 0, and then letting the held be constant, one gets the 
how equations for all n-point functions in a constant background held 0. These equations 
can be represented diagramaticaly by one loop diagrams with dressed vertices and propa¬ 
gators (see e.g. [19]). For instance, the how of the 2-point function in a constant external 
held reads: 


-p; </>) = 




xG'^((? + pf-, + R -R 0) 

-P, g, -g; 0)G^(g^ 0 )| , 

( 5 ) 

where 


G'«^(g^0) = F®(g,-g;0) + R^(g2). 

(6) 

The corresponding diagrams contributing to the how are shown in Fig. 1. 



Flow equations for the n-point functions do not close: for example, in order to solve 
eq. (5) one needs the 3- and the 4-point functions, F(f) and F[)^^ respectively. At this point 
we observe that because of the shape of the regulator, only internal momenta q smaller 
than K contribute to the how, i.e., to the integral in the r.h.s. of eq. (5). One refers to this 
property as to the decoupling of high momentum modes. Besides, the regulator insures 
also that all vertices are smooth functions of momenta. Suppose then that one wants to 
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calculate n-point functions at small external momenta, say, pf ^ k^. Then all momenta 
are small. This, together with the fact that the n-point functions are smooth functions of 
the momenta, make it possible to expand the n-point functions in the r.h.s. of the 
flow equations in terms of and or equivalently in terms of the derivatives of 

the held. 

Such considerations are at the basis of the derivative expansion, which is usually 
formulated in terms of an ansatz for the effective action, rather than in terms of an 
approximation for the n-point functions. Its zero order, the LPA, assumes that the effective 
action has the form: 

= J d^x + Vfc(0)|. (7) 

The derivative term here is simply the one appearing in the classical action, and 14(0) is 
the effective potential. The how equation for 14 is a closed equation that is easily obtained 
by assuming that the held 0 is constant in eq. (3): 

<9«I4(0) = ^j 0), (8) 

with G^^^{q‘^;(f)) given by Eq. (6) in which 

0) = + <9V«/(90^ (9) 


as obtained from Eq. (7) 

Higher order corrections to the LPA include terms in the ehective action with an 
increasing number of derivatives. Although there is no formal proof of convergence, the 
expansion exhibits quick apparent convergence if the regulator Ri^{q^) is appropriately 
chosen [17,20, 21]. In practice, the LPA reproduces well the physical quantities dominated 
by small momenta (such as the ehective potential or critical exponents) in all theories 
where it has been tested (see, for example, [5,7]). Higher order corrections improve the 
results. The expansion has been pushed up to third order [21], yielding critical exponents 
in the Ising universality class as good as those obtained with the best accepted methods. 


3 A new approximation scheme 


The derivative expansion, as described above, strictly makes sense only for momenta 
not much larger than \JiR + m'f., where is the running mass. Thus, at criticality, and 
in the physical limit k —0, it provides information only on n-point functions and their 
derivatives at zero momenta. However, by focusing on the n-point functions rather than 
the effective action, we can generalize slightly the arguments on which the derivative 
expansion is based in order to set up a much more powerful approximation scheme. We 
observe that: i) the momentum q circulating in the loop integral of a flow equation is 
limited by k] ii) the smoothness of the n-point functions allows us to make an expansion 
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in powers of independently of the value of the external momenta p. Now, a typical 

n-point function entering a flow equation is of the from + q,Pn — Q] 0), 

where q is the loop momentum. The proposed approximation scheme, in its leading order, 
will then consist in neglecting the g-dependence of such vertex functions: 

rJf)(pi,P2, -.Pn-l + q,Pn- g; <t>) ~ rjf)(pi,p2, ...,Pn-l,Pn, 0 ). (10) 

Note that this approximation is a priori well justihed. Indeed, when all the external 
momenta Pi = 0, eq. (10) is the basis of the LPA which, as stated above, is a good ap¬ 
proximation. When the external momenta pi begin to grow, the approximation in eq. (10) 
becomes better and better, and it is trivial when all momenta are much larger than k. 
With this approximation, eq. (5) for instance becomes: 


Sk^kKp, -p] <(>) = / 0)r(fnp, 0, -p- 0) 

xG^((g + pf- 0)r® (-p,p, 0; 0)G^(g^; 0) 

-]^G^{q^-(t))T^^\p, -p, 0, 0; (t))G^{q^- 0)| (11) 

Note that we do not also assume g = 0 in the propagators. The reason for this will become 
clear shortly. 


Now comes the second ingredient of the approximation scheme, which exploits the 
advantage of working with a non vanishing background held: vertices evaluated at zero 
external momenta can be written as derivatives of vertex functions with a smaller number 
of legs (this observation can be found, in a different context, in [1] and in [22]): 


rir^^nPi,P2,--,Pn,O;0) = 


<9rl"^(pi,P2,--p«;0) 

dcj) 


To prove this relation, expand around an arbitrary constant held 0°: 


( 12 ) 


/ d‘‘x,...d'‘xnl4>{xi) - - ’t’°]rp{Xu--,Xn;4’o)‘ (13) 

n n. J 

Since T^f^j does not depend on 0°, dV Jycl)]/dcff' = 0. Taking the derivative of eq. (13) with 
respect to 0o one then gets: 




dTG)(xi,X2, ...Xn,0° 
( 90 ° 


(14) 


from which eq. (12) follows after a Fourier transform. By exploiting eq. (12), one easily 
transforms eq. (11) into a elosed equation (recall that and T^.^^ are related by eq. (6)): 
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Fig. 2. The two diagrams contributing to the flow of the n-point vertex and containing vertices 
with more than n legs. 
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(15) 


It is interesting to emphasize the similarity of this equation with eq. (8): both are 
closed equations because the vertices appearing in the r.h.s. have been expressed as deriva¬ 
tives of the function in the l.h.s.. This is the main result of this paper. 

The approximation scheme presented here is similar to that used in [12]. There also 
the momentum dependence of the vertices was neglected in the leading order. However 
further approximations were needed in order to close the hierarchy. The progress realized 
here is to bypass these extra approximations by working in a constant background held. 

The construction of closed equations for the n-point functions with arbitrary n (n > 3) 
follows the same lines as that of the equation for the 2-ppoint function. The how equation 
for involves all m-point functions with m < n+2. Since the r.h.s. of the how equations 
have the structure of one loop integrals, the contributions involving and are 

of the types shown in Fig. 2. When the loop momentum entering these vertices is taken 
to be zero, in line with eq. (10), the vertex of order n -|- 1 has one leg at zero momentum 
and the vertex of order n -|- 2 has two; thus, according to eq. (12), they can be written 
as derivatives of the n-point function It follows that the equation for T)”^ is a closed 
equation, assuming of course that all other functions T)™^ with m < n are determined 
similarly. 

One can go beyond the leading order approximation, based on eq. (10), in the following 
way. Focusing on eq. (5) for F[2)(p, —p; 0) one can solve simultaneously the how equations 
of F)20 F(f^ and F](^0 with no approximation in the how equation for F(20 but using (10) 
in the right-hand-side of the how equations for F|.^) and F(fb In this way one would de¬ 
termine F).^) and with a “leading order” precision and F(2) with a “next-to-leading 
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order” one. By iterating the procedure, which amounts to including more equations, one 
would get better approximations for a larger number of n-point functions. It is easily 
verihed that such an iterative approximation scheme reproduces perturbation theory at 
high momenta: for instance, if the bare vertices of the theory are momentum independent, 
the leading order approximation contains the exact one loop (this is only true if we keep 
q in the propagators, as we did in (11)). More generally, a simple analysis shows that to 
get the expression of a given n-point functions at order m loop, one has to consider all 
equations up to that for the n -|- 2m — 2-point function. (In a theory with momentum 
dependent bare vertices, perturbation theory is also reproduced, but more flow equations 
are required to match perturbation theory at a given number of loops.) Of course, such a 
scheme is not only accurate at high momenta, where it reproduces perturbation theory, 
but also in the small momentum, possibly critical, regime, where its accuracy is at least 
comparable to that of the derivative expansion. In practice, the scheme that we have de¬ 
scribed may become rapidly cumbersome. However, it can be simplified by using a further 
approximation consisting in a systematic expansion in of the n-point functions in 

the right-hand-side of flow equations. This further approximation preserves the essential 
property of leading to a closed set of equations. We have solved the flow equation for 
the 2-point function in leading order and obtained results which are comparable to those 
obtained in [12]. More complete numerical studies will be presented in a forthcoming 
publication. 


4 The 0{N) model in the large N limit 


We turn now to the 0{N) model in the large N limit, where analytical results can 
be obtained. Our goal here is twofold. First, to prove that, in its leading order, our 
approximation is in fact exact in the large N limit. Second, to solve analytically the 
NPRG equations and hnd the n-point functions in a direct way, i.e., without introducing 
an auxiliary held, as commonly done [23]. 

The large N limit of the 0{N) model within the NPRG has been thoroughly studied 
in [24], where it is shown that the effective action takes the form: 

p.[0] = I d<^x + f.[p], (16) 

where p{x) = (f)°‘{x)(f)°‘{x)/2, a = 1,... ,N. Using this expression, the authors of Ref. [24] 
show that the LPA how equation for the ehective potential is exact. Moreover, they hnd 
an interesting analytical solution of the equation. However, the use of (16) in the how 
equation of a n-point function at zero held does not give a closed equation [25]. In this 
section we obtain a system of closed equations that can be solved analytically. 

We shall again focus, for the sake of illustration, on the 2-point function. Its exact 
how equation in an external held is the generalization of eq. (5)) to the 0{N) model: 



( 17 ) 


^n'^abiP: -p;0) = / {Gijiq'^; 0)r3i(p, g, -p - g;«; 0 ) 

xGkiiiq + p)^;«; 0)rwi(-p>p + g^ -g; 0)G'm*(g^; 0 ) 
-ijGij{q^- K- 0)ritifc(p, -p, g, -g; (fyGkM^- k; 0)|. 


At this point, it is useful to recall the structure of the n-point functions of the 0{N) 
model in the large N limit. These are obtained by taking the functional derivatives of 
(16), letting the background held to be constant, and then taking the Fourier transform. 
We obtain in this way the following expressions for the 2-, 3- and 4-point vertices: 


rl?(p, -p; K; 0) = (p^ + V^{p))5ab + 0a0bfi^nP, -p; p), (18) 

rl?c(PnP2,P3; -Pi;p) + 065acfi^HP2, -P2;p) 

+ 0Abf(,^nP3, -P3;p) + 0a060cf)fHPnP2,P3;P), (19) 

rilL(Pl>P2,P3,P4;K;0) = +P2,-Pl - P2] p) 

+5ac5bd^^i^\Pl + P3, -Pl - P 3 ; P) + ® (Pl P 4 , “Pi “ P 4 ; p) 

+^afe0c0df®(pi -hp2,P3,P4;p) + ^ac060(if ® (pi Ps, P2, P 4 ; p) 

+5ad060cf®(pi -hp4,P2,P3;p) + ^fec0a0df ® (p2 + P3, Pl, P4; p) 

+56d0a0cf®(p2 +P4,P1,P3;P) + ® (pS P 4 , Pl, P2; p) 

+0a0fe0c0df^'^HPl,P2,P3,P4;p) • (20) 

As for the propagator, it can be written in terms of its longitudinal and transverse com¬ 
ponents: 

G'ab(p^;0) = G't(p^;p) +Gl(p^;p)^|^, (21) 

with 


Gt{p^\k; p) =p^ + V^{p) + Rn{p^) , 

G'Z^(P^; P) =P^ + K(p) + 2pf(2)(p, -p; p) + i?^(p2). (22) 

We can now replace these expressions in eq. (17) and keep only the leading terms in 
1/N. To do so, notice that the non trivial large N limit of the 0{N) model is obtained 
when the effective action T^ is considered to be of order N and the 0 held of order \/N 
(see for example [24]). Then, the 2-point function is of order 1, the 3-point function of 
order 1/V^ and the 4-point function of order 1/N. Simple counting of N factors in the 
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flow equation (17) gives a right hand side of order 1/N . Thus, the only way of having 
both sides of order 1 is to collect an explicit N factor from the trace of an identity matrix. 
The surviving terms yield; 

-p; K]<i))=N J 

+<Pa<Pb -P: p)y Gt{{p + qf] K- p) - ® (p, -p, 0; p)| . (23) 

Now comes the important observation: the vertices contributing to the flow in the 
large-iV limit are g-independent. That means that our leading order approximation be¬ 
comes exact in the large-iV limit of the 0{N) model, as anticipated. Notice that for this 
to be true, we need to keep the q dependence of the propagators in the flow equations. 
However, in the large-N limit, the momentum dependence of the transverse propagator 
(the only one appearing in (23)) is simply the bare one, as shown in eq. (22). 

Substituting in eq. (23) the expression (18) of the 2-point function and performing 
the tensor decomposition, one obtains the system of equations: 

At r rl^n 

'9«K'(p) = -yK'(p) J k] p) 

-p; p)=N J -^^d^R^{q‘^)G^q‘^; k; p) 

X {(ff (P. -p; ^Griip + p) - 1 ^T>(P,-P;P) j p,, 

where we have used, as in eq. (12), 

ff(0,0;p) = K:'(p). ff(p,-p,0;p) = ^!2^k^=^, (25) 

which allowed us to close the equations. In the equations above, V^{p) and V”{p) denote 
respectively the hrst and second derivative of the effective potential with respect to p. 

The hrst of eqs. (24) is simply the derivative of the equation for the potential in 
the large-N limit which was solved analytically using two different methods [24,26]. The 
second equation, which, to our knowledge, has not been presented before in the literature, 
can be solved also using any of these two methods. Here we present the solution using the 
method of [26]. 

To this aim, we set W = V]((p)) use (k, hh) as independent variables instead of 
(k, p). In order to do so, we introduce the inverse function relating p to IT: 

p = ( 26 ) 
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so that: 


dJ.iW) = -f^{W)d.W{p), f'M) = l/W'{p), (27) 

where /' and W are the derivatives of / and W with respect to their respective arguments. 
Making the change of variables in eq. (24), and using the relations (27), one gets: 





N f 




-p]p 


-1 


= -N 


{27rY {q^+ W + R^{q^)y 

f d'^q d^R^{q^) 


{2TrY (g2 + pp + Rniq'^W (q + p)'^ + w + R^{{q + p)^) 

(28) 


At this point we observe that the right hand sides of both the expressions above are total 
derivatives: 


r d‘^q d^R^{q^) ^ f d^^q 1 

J {2ttY {q^ + W + i?«(g2))2 “ (2ttY q^ + W + R^{q^) 
r d'^q d^R^jq^) _ 1 _ 

J {27rY {q^ + W + i?.(g2))2 (g + p)2 + p/ + + p)2) 

f _ 1 _ 1 _ 

2 " y (27r)'' q^ + W + R^{q^) {q + py + W + R^{{q + p)^)' 

(29) 

This allows us to integrate analytically the equations. For the potential, one obtains the 
well known result [5,26]. By imposing the initial condition that for k = A the 4-point 
function is the bare vertex, we get from eq. (20) 

ri=A(p,-p;p) = (30) 

corresponding to a 0^ interaction in the classical action, of the form {u/4\){(j)a4>a)‘^- We 
then obtain the longitudinal part of the 2-point function (in the limit A —cx)): 

(P,-P;P) = T . I Vn r 1 ^ 1 (^1) 

-‘■■^6 J q^ + W+R^{q‘^) (q+pp+W+R4{q+pp) 

This is the usual sum of bubbles in large N, which can be found using other methods [23]. 
Notice however that our solution follows directly from the NPRG without introducing 
any auxiliary held. The dependence on W rehects the presence of the background held: 
W plays here the role of a (held dependent) “mass term”. 

The same method can be applied to the other n-point functions. The calculations are 
straightforward, but the explicit results involve lengthy expressions. Let us just give here 
the result for the 3-point function: 
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fL'Hpi,P2,P3; 1^) = ivf(2)(pi, -p,- W)t^^\p2, -P2; W)t<;^\ps, -P3-, 


X 


d^q 

{2'kY 


p)GT{{q + Pi)^; «; p)GT{{q — Pzf] i^', p)- 


( 32 ) 


5 Conclusions and perspectives 


In this paper we have presented an approximation scheme to solve the NPRG equa¬ 
tions and obtain the n-point functions for any external momenta. Our proposal is as gen¬ 
eral and systematic as the derivative expansion: results can be systematically improved 
as one can write, at any order, a closed set of equations. We have provided here, as an 
illustration, an application to the 0{N) model in the large N limit: there, the method 
turns out to be exact at leading order, and it provides an economical procedure to obtain 
the analytical expressions of the n-point functions at arbitrary momenta. Clearly, how¬ 
ever, the method is more general. It is based on an approximation that can be exported to 
other held theories. Work on such extensions, in particular to gauge theories, is underway. 
This, as well as numerical studies for physically interesting problems, will be reported in 
forthcoming publications. 
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